Scaling of transition temerature and Cu02 plane buckling in the cuprate 
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A universal feature of all high temperature superconductors is the existence of a chemical com- 
position that gives maximum Tc, separating so called under and over-doped region. Recently it is 
seenci that both Tc and the buckling angle of the Cu02 planes goes through a maximum at the same 
doping level. We show that only for optimal doping concentration the Fermi surface touches the 
M(0, tt) point in the BZ, where the matrix element for interlayer pair tunneling amplitude is largest, 
so that the gain in delocalization energy by tunneling(in pairs) along the c axis is largest. Buckling 
of the planes on the other hand modulates the separation between the planes and thereby modulates 
the interlayer pair tunneling amplitude. That is why both Tc and buckling angle( Oxygen atom dis- 
placement out of the plane) scales the same way with doping concentration. We have calculated Tc 
and buckling angle for various doping concentration. The agreement with experiment is remarkably 
good. We also point out the possible reason for large(about 1 percent) change of the buckling mode 
phonon frequency, accross the tj*jansition temperature for optimal doping concentration, observed in 
neutron scattering experimental 

PACS : 63.20.Kr, 74.20.Mn,74.25.Kc,61.12 Ex, 4.72.Bk 



I. INTRODUCTION 



Since the discovery of High Tc cuprates a wide variety of experiments have shown strong evidence of a link between 
anomalies in lattice structure and the onset of superconductivity in these materials. For a review of experiments 
and analysis please seeEI. Ion channeling experimentsu show anomalous changes in the oscillation amplitudes of Cu 
' and O atoms at and above Tc. Ultrasonic measurementaJ shows anomalies in elastic constants. Infrared reflectivity 
measurementfu of high quality samples have shown clear connection between c-axis phonon modes and the charge 
I . carriers. _ 

' Recently it is seenEl that both Tc and the buckling angle of the Cu02 planes goes through a maximum as the doping 
Ch , is varied. The buckling corresponds to a static displacement of the oxygen atoms out of the plane, with Ui±x = —Ui±y, 
O ■ where Ui±x is the displacement of the oxygen atoms at i ± x along the c axis, i being the Cu atom position. 

In inelastic neutron scsaterring experimentsB it is also seen that the buckling mode phonon at 340 cm ~^ ( called 
Aig ) softens by more than 1 % accross the transition temperature. On the other hand, other modes having frequency 
close to the 340 cm^^ mode like ( 302 cm^^ ) and E2U ( 347 cm^^) and Big (285 cm are studied in detail by 
inelastic neutron scattering experimentg3. The Aig mode corresponds to the vibrations of the planar oxygen atoms 
out of the plane, and out of phase with the adjacent CuO plane (buckling mode). The Big ( here the oxygen atoms in 
two adjacent planes vibrate in phase and out of the plane), and the i?2u mode£ here the atoms vibrate in the plane ) 
shows much less softening ( estimated to be less than .4% and .3% respectivelyo) . This compares well with the optical 
data though for the £'211 there is no optical data available to compare with. 

In this paper we shall be exploring the phenomena of static buckling of the planes as well as,the the phonon mode 
softening within the interlayer tunneUpg mechanism of Wheatley, Hsu and Anderson(WHA)Ll and its more refined 
version proposed by Chakraborty et alH. We will show that Tc and the buckling angle varies with doping concentration 
the same way. This and the large softening of the buckling mode phonon can all be understood in a unifying way. 
The main ingradient of this mechanism is that the normal state is a non Fermi liquid and hence the single particle 
hopping out of planes becomes incoherent. However coherent propagation of pairs of particles becomes favourable 
at low enough temperatures leading to coherent pair fluctuation and the consequent superconductivity. The induced 
Josephson coupling between the planes, or the amplitude for pair tuneling is quadratic in (bare single particle 
hopping amplitude between the planes). 

Any static displacement or dynamic fluctuations of the atoms in the unit cell that modulates the distance between 
the planes effectively modulates the pair tunneling amplitude between the planes. The main thesis of our paper, 
being, since the Aig mode effectively brings the planes closer together leading to a decrease of the hopping paths for 
the electrons between the planes. This in turn makes the Josephson coupling stronger, the electronic system gains 
energy by condensation and the corresponding optical phonon mode(g = 0) frequency gets softer, leading also to an 
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average static buckling. We find that the greater part of the large frequency shift observed for the Aig mode(buckling 
mode) is due to the screening by interlayer pair propagation , or due pair current-phonon coupling. The Big mode 
on the other hand, being in phase oscillations of the atoms keeps the distance between the planes unchanged and so 
it does not couple with the pair-current flow accross the planes. 

II. INTERLAYER TUNNELING MODEL 

We begin by writing the gap equation from interlayer tunneling^l, 

A. = T.jik)^t^nh^ + ± V,,, ^tanh^, (1) 

This equation can be obtained by considering two close Cu-0 layers as in BI 2212 coupled by a Josephson tunneling 
term of the form 

k 

where t_i{k) is the bare single electron hopping term between the two coupled layers c and d and i is a band structure 
parameter in the dispersion of electrons along the Cu-0 plane. Finally, Tj{k) in the right hand side of equation (1) 

is given by Tj{k) = ^^^^^ . where tj_(fc) = ^(cos(A:a::) — cos(fcj/))^. The dispersion of electrons along the Cu-0 plane 
is chosen to be of the form 

t{k) ~ — 2t ( coskx + cosky ) + At' coskx cosky — 2t" {cos2kx + cos2ky), (2) 

with t = 0.25 eV ,y = 0.45 , ^ = 0.2. We also choose ep = -0.45 eV correspDnding to a Fermi surface which is closed 
around the T point. These choices are inspired by band structure calculationstl. Note that the Josephson coupling term 
in Hj conserves the individual momenta of the electrons that get paired by hopping across the coupled layers. This 
is as opposed to a BCS scattering term which would only conserve the center of mass momenta of the pairs. This is 
the origin of all features that are unique to the interlayer tunneling mechanism. This term has a local U{1) invariance 
in fc-space and cannot by itself give a finite Tc. It needs an additional BCS type non local interaction in the planes 
which could be induced by phonons or residual correlations. Here we assume the inplane pairing interaction to be 
d-wave kind i.e, Vkk' = Vq fkfk' with fk — cosk^ — cosky. Tj{k) can be infcrcd from electronic structurecalculations. 

As shown in reference!, it is adequate to choose t±{k) = -^(cos k^ — cos fcy)^. According to Anderson, it is the fc-space 
locality that leads to a scale of Tc that is linear in interlayer pair tunneling matrix element. He finds that in the limit 
Tj > Vkk', ksTc ^ and in the other limit, fesT^ « hujoe~po^ , where too Po and Vq are Debye frequency, density 
of states at the fermi energy, and fermi surface average matrix element Vkk' ■ The important point being that even 
with a little help from the in plane pairing interaction the interlayer tunneling term can provide a large scale of Tc. 
In our naive analysis with a fixed cut off ljd we find that the Tc drops when the fermi surface moves away from the 
M(0,7r) point in a parabolic fashion. This because of (1). the density of states is large near the M points, and only 
for the optical concentrations, the fermi surface touches these regions, and (2) The momentum dependence of both 
the inplane BCS pairing interaction(d-wave kind) and the interlayer pair tunneling amplitude are such that they have 
large values near the M also. So in a way the Tc going through a maximum near the optimal doping concentration, 
is expected for both the interlayer tunneling mechanism and as well as for ordinary d-wave superconductor. Now we 
shall discuss the variation of the buckling angle with the doping concentration. 

Some other feature of the WHA mechanism are, (1). The gap values near the F — X(0, ±7r) points are very much 
temperature insentive at low temperatures, and drops rapidly to zero near the Tc. It is also very sensitive to small 
doping variations, near the optimal doping concentrations, i.e as the fermi surface moves away from the T — X points 
because it cannot then take advantage of the large pair tunneling amplitudes near those points in the BZ. (2). The 
gaps along the diagonal direction in the BZ(F — M directions) are on the other hand very sensitive to temperature 
variations, falling faster than BCS with temperature. It is also insensitive to small doping variations near the optimal 
doping concentrations. We shall discuss the possible consequences of these features of the WHA gap later. 
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III. BUCKLING OF COPPER OXYGEN PLANES 



The interlayer pair hopping term for a bilayer system is 

+•2 



E^_l(^) / It It 2 2 ,1, \ 



k 

Now the interlayer hopping matrix element t±_ will be modulated due to the vibration out of plane of the oxygen 
atoms. We emphasize here again, that for the Big mode the atoms all vibrate out of the CuO planes but the vibrations 

in two adjacent layers arc in phase and hence the effective tunneling matrix clement of the Zhang-Rice orbitals (a 
linear combination of the oxygen p orbitals) from plane to plane dont change, but for the Aig mode t±^ will change. 
We can do a Taylor expansion of t±_ in small oxygen atom displacements. ti_(u) = t±{Q) + a{^)u where c is the z 
direction layer separation and a is a constant of the order of one. Putting this back in the interlayer term we get the 
extra piece, 

U"" y^ t\{k) u It 2 „2 , /, „ ^ , ^ u It 2 „2 , , „x 

k k 

With these modifications the total effective Hamiltonian for the phonon modes corresponding to oxygen vibration out 
of plane, will be given by: 



where we have taken the anomolous averages. As we can see that the second term renormalises the frequency of the 
buckling mode phonon, and the term linear u shifts the harmonic oscillators without affecting its energy(this will 
correspond to static displacement of the oxygen atoms causing buckling of the Cu02 planes). 

The renormalised frequency of the buckling mode phonon wofe and the static displacement of the oxygen atoms out 
of the plane uq is given by, 

^, , 2A 2acA 

and Uo = ——2 7- • (4) 



where, A = {{c%c\) {c\,c%) + h.c.) 

We take w° to be equal to 340 cm ^ or 10.2 THz for fc = 0. In Fig: 1 and 2, we show the doping variation of both 
the superconducting Tc and the buckling angle 6 = tan~(2uo/a) where a is the latice spacing(Cu-0-Cu distance). We 
see that both Tc and angle 6 scales the same way with the doping concentration. As we have seen before, due to the 
momentum dependence of the pair tunneling amplitude at the optimal doping concentration, large number of pairs 
can delocalize by tunneling from plane to planes (thereby reducing the kinetic energy). The buckling mode phonon 
helps in this process. That is why the buckling mode phonon frequency softens and also develops a static buckling 
distortion. 



IV. PHONON MODE SOFTENING 



Now let us take a look at the phonon frequency softening below Tc for different phonon modes . We shall specifically 
look at the Aig and Big modes( corresponding to out of phase and inphase vibrations of the oxygen atoms at adjacent 
bilayers). 

Superconductivity induced change in frequency and linewidths have been used to estimate the superconducting 
gap function. This is accomplished by determining the phonon self energy changes that occur below the Tc- The 
imaginary part of the phonon self energy is proportional to the phonon linewidth (inverse lifetime) which is measured 
as a function of temperature. If the phonon energy is less than twice the gap energy , it cannot decay by pair breaking 
and thus a decrease in linewidth below can result because of a decrease in the number of available decay channels. 
If the phonon energy is greater than the gap, it has sufficient energy to dissociate an electron pair and the number of 
possible decay routes can increase below Tc, resulting in an increase in linewidth. 

The real part of the self energy is proportional to the phonon energy and the superconductivity induced changes are 
determined from a measurement of the phonon frequency as a function of temperature below Tc- Qualitatively one 
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expects that, for phonon energies below or near the gap , the electron phonon interaction will push the frequencies to 
smaller values. Conversely for larger energy phonons (w > 2A) are expected to undergo a hardening. 
The phonon self energy, assuming a polarisation bubble with no vertex correction is given by, 

T 

E('^' *^") = ]^ E E l^7(^, k + g)pTr(r3G[fc + 

7 km 

where ra is a pauli matrix, G{k, iujn) is the fully interacting green's function, uj„ ii^n) is the fermion(boson) Matsubara 
frequencies, and A^(fc, k + q) is the electron phonon matrix element for scattering an electron of momentum k to k + q 
with momentum transfer q to or from a phonon with branch index 7. The standard electron phonon interaction being, 

He-ph = E KilWk-q.a'^k,aialj + a^q,^) (6) 
k,q,-y.(T 



- efcT3 + AfcTi 



In BCS theory the superconducting Green's function is, 

G{k,iuJn) ^ 

Substituting this in the earlier equation and taking g = we get, 

where Ek = (e| + A|)^/^ and X^k = X-yik, k). Raman shift is the difference between the superconducting and normal 
state self energies, where normal state self energy can be extraced from the above equaition by putting = 0. One 
still has to analytically continue ( ivn —fV + iS and the real part describes the frequency shift and the imaginary part 
describes the change in phonon linewidth upon entering the superconducting state. 

The simplest method of evaluating it, is to numerically evaluate it on the lattice after solving the self consistent gap 
equation. We have done so iq this paper. We borrow the estimates of the electron phonon coupling function from the 
estimates by Normand et alliJ, i.e. t{u) — t{0)[l — 2.03^] where t is the inplane hopping matrix element, u is the small 
oscillation amplitude of the oxygen atoms out of the plane, and a is the lattice parameter. Also for the Aig mode 
symmetry (uf = —u^ ) and for the Big mode { uf = ). The electron phonon coupling function Afc^^ for Aig and 

phonons will be Xkk — 4.06-^^^ 2)mU^(''')' '^^^'^re j{k) = cos fc^ ± cos /cj^ for Aig and Big phonon mode. The constants 
m and to are the mass of the oxygen atom and the bare phonon frequencies((7 = 0) in the normal state for the Aig 
and i3ig(340 cm~^ and 307 cm~^) modes. This contribution to phonon frequency softening is due to the screening 
by the charge carriers in the planes. We have already calculated the phonon mode softening due to interlayer pair 
hopping process (additional screening due to pair tunneling process). 

Both the contributions that shifts the phonon frequencis, the usual self energy contribution and the screening effect 
due to the pair propagation along the c-axis are shown as a function of temperature for the Aig and Big phonon 
modes in Fig. 3 and 4. 

For the Big phonon the second contribution will be missing and is plotted in Fig. 4. We see that the pair tunneling 
contribution is much larger than that of screening by inplane charge carriers. Seecndly from the form of electron 
phonon coupling function A^fc as well as the pair hopping contribution(the k dependence of t±{k)), it is clear that, the 
contributions to the frequency renormalisation of Aig and Big phonons comes mostly from the M(0,7r) points, and 
from the points along the diagonal respectively. This in the light of the WHA gap nature would mean that the change 
in frequency accross the Tc will be more abrupt for the Aig phonon, while Big frequency change will be gradual. At 
lower temperatures Aig frequency willi very sensitive to small doping variations near the optimal filling, decreasing 
as the fermi surface moves away from 

In conclusion, we find that, both scaling of Tc and the buckling angle with doping variation and large frequency 
softening of the buckling mode phonon compared to other optical modes of comparable frequencies, can be understood 
with the WHA mechanism. i_. 

In La2-x-ySrxNdyCuOii compounds, it is noticedEJ that with increase of Nd concentration, the buckling(or the 
tilting of Cu)q octahedra) increases while the Tc falls, i-c for the same doping concentration, the Tc falls with buckling 
angle. It was suggested by Bonesteel, Rice and Zhangt3 that the charge carriers gets coupled with the buckling angle 
via spin orbit coupling. It was shown by Bonesteelllj that spin orbit scattering can be a strong pair breaker provided 
tye material is in a structurallydisordered phase, in which the buckling distrotion is coherent on small length scale but 
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is completely disordcd on large length scales( i.e. a spatial variation of buckling angle unlike what we have considered). 
It is possible that this series of materials have a lot of unrelaxed structural disorder. 
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Fig.l Transition temperature versus doping concentration in the interlayer tunneling model. 
Fig.2 Transition temperature versus doping concentration in the interlayer tunneling model. 

Fig. 3 Frequency versus Temperature of Aig phonon for three different chemical potentials(/x = -457.0,-450.0,-440.0 
meV, i.e doping concentrations of 0.19, 0.18, 0.16). The corresponding T^s are 92.5, 88.0 and 85.0 degrees. The 

upper three curves shows the screening of inplane carriers (the; bare polarisation bubble) , the middle ones are the 
contributions due to interlayer tunneling, and the bottom ones are the total contributions. 

Fig. 4 Frequency versus Temperature of Big phonon for three different chemical potentials(/i = -457.0,-450.0,-440.0 
meV, i.e doping concentrations of 0.19, 0.18, 0.16). 
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